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1. Bk 

The size of indirect CP violation in the neutral Kaon system is, in experiment, parameterized 
by £, which is proportional to the kaon bag parameter Bk defined as 

^ ^ {K^m^{^-y^)d][sy^{i-y5)d]\Ko) 

|(^ok"7M75^/|0)(Ok-7M75^/|^o) 

A precise determination of constrains the CKM matrix, which might lead us to a window of 
new physics beyond the standard model. Hence, there have been a number of efforts to calculate 
Bfi with higher precision. As pointed out in Ref. [|]], the large scaling violation observed in the 
calculation using unimproved staggered fermions can be reduced remarkably by improving stag- 
gered fermions with HYP fat links. Even though the scaling violation is taken care of by improving 
staggered fermions, there has been an uncertainty originating from the quenched approximation. 
Hence, it has been crucial to perform a numerical study in unquenched QCD so that we can re- 
move the uncertainty due to quenched approximation. In this paper, we describe our first attempt 
to perform a numerical study on B^. in partially quenched QCD, while minimizing the scaling vio- 
lation using improved staggered fermions. We use HYP staggered fermions as valence quarks and 
AsqTad staggered fermions as sea quarks (we call this a "mixed action"). Details of the simulation 
parameters are summarized in Table |l|. 



parameter 


value 


sea quark 


AsqTad Staggered 


valence quark 


HYP Staggered 


i8 


6.76 {Nf = 2+1 QCD) 


# of confs 


593 


lattice 


20^ X 64 


sea quark masses 


^u4 = 0.01, m-s = 0.05 


valence quark mass 


0.01,0.02, 0.03, 0.04, 0.05 



Table 1: Parameters for the numerical study 

We measure weak matrix elements and hadron spectrum over a subset of the MILC gauge 
configurations |^. We set up U(l) noise sources at T = and T = 26, which project out only 
pseudo-Goldstone pions (75 (g) ^5) and exclude all the other non-Goldstone pions. The kaon signal 
as a function of Euclidean time T is shown in the lefthand side of Fig. [l| Here, we observe that 
there is a noticeable contamination from excited states in the range of < T < 10. Hence, the best 
fitting range for Bk islO<f<15in order to exclude the contamination from excited states, as you 
can see in the righthand side of Fig. 

Ref. [^] presents a chiral behavior of Bk and its finite volume effect in the case of Nf = 
2 partially quenched QCD. Although this result is interesting, it does not apply directly to our 
numerical study mainly because it is in Nf = 2 + I partially quenched QCD. Recently, Van de 
Water and Sharpe have calculated the chiral behavior of Bk in Nf = 2+1 partially quenched QCD 
using staggered chiral perturbation theory [Q]. The result for the degenerate valence quark mass 
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Figure 1: niK vs. T (left) and Bk vs. T (right) at —md= 0.03 



combination {m^ = niy) is 



1 



48712/^ 

2 

xy 



^conn ~l~ ^disc 



6mU{ml)-m{m' 



xy) 



^disc = 

/(X) = Xlog(X/A2)+F.V. 
r(X) = -[Iog(X/A2) + l]+F.V. 



(1.2) 

(1.3) 
(1.4) 
(1.5) 
(1.6) 



where / = 132 MeV, c,- are unknown dimensionless low-energy constants and F.V. represents a 
finite volume correction. The notations for the various meson masses are as follows for those 
composed of sea quarks: 



and as follows for those composed of valence quarks 



my 



2/lm 



m 



■xy 



■ n{mx + my) 



(1.7) 



(1.8) 



Here, we set sea quark masses to m„ = 7^ m^. and the two valence quark masses are nix and my. 
In this paper, we consider only the case of ntx = my {i.e. m\ =mY = m^). 

In Fig. |2[ we plot Bf^ data as a function of kaon mass squared M\. Here, all of them have 
degenerate valence quarks {mx = my). We fit the data to the form of Eq. (L2) suggested by chiral 
perturbation theory: 



ci 1 + 



1 



48712/^ 



+ C2mly + c^mly 



(1.9) 
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Figure 2: Bk vs. M| 

where the cut-off scale A is set to A = Anf and the remaining scale dependence is absorbed into 
C2. The fitting results are summarized into Table ^. 

In Table ^ the is rather high when we take into account the fact that all the data are 
correlated. In fact, one can observe that the fitting curve does not fit the lightest two data points 
very well in Fig. ^. In other words, the fitting curve miss the lightest two data points in the opposite 
direction. In Ref. [Q], it is pointed out that the contribution from the non-Goldstone pions are so 
significant that the curvature of the fitting curve becomes smoother, which is consistent with what 
we observe in Fig. ^ However, the full prediction from staggered chiral perturbation contains 21 
unknown low-energy constants for a single lattice spacing and 37 unknown low-energy constants 
for the full analysis In order to determine all of them, it is required to carry out a significantly 
more extensive numerical work including data with non-degenerate quarks 



parameters 


average 


error 


Cl 


0.4488 


0.0162 


C2 


-1.4883 


0.1883 


C3 






C4 


1.4533 


0.2207 


X^/d.o.f. 


0.6812 


0.6 



Table 2: Fitting results for Bk 

In the current analysis of Bk data, we match the lattice results to the continuum values at 
the tree level. In this respect, the results are preliminary. Hence, we plan to calculate the one-loop 
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corrections to the four fermion operators in near future. In addition, note that there exists a different 
approach to Bk using staggered fermions [§]. 



2. 5(3/^) and 5(3/2) 

In the case of direct CP violation e'/e, there are two contributions which interfere with each 
other destructively: the A/ = 1/2 part and A/ = 3/2 part. The A/ = 3/2 part is dominated by the 
electroweak penguin contribution from Q$ and Qj. In this section, we focus on the bag parameters 
of Qs and Qj. 
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Figure 3: Bj vs. T at m,. = = 0.03 (left) and Bj vs. (right) 



,(3/2) 



The B) ' ' data as a function of Euclidean time T at quark mass irix = 0.03 is presented in the 
lefthand side of Fig. |3[ The fitting range is the same as for Bk- In the righthand side of Fig. ^, B^j^^"^ 
data is plotted as a function of M|-. We fit the data to the form suggested by chiral perturbation 
theory: 



B 



(3/2) 



ci + C2mly + c^rn^ log(m^./A- 



■xy- 



(2.1) 



The results are summarized in Table ^. Here, note that we match the lattice results to the continuum 
values at the tree level. In this respect, the results are preliminary. 



parameters 


^(3/2) 


r(3/2) 


Cl 


1.3068 ± 0.0078 


1.2505 ±0.0074 


C2 


-0.1889 ± 0.0093 


-0.1439 ±0.0091 




0.4080 ± 0.0228 


0.3175 ±0.0222 



Table 3: Fitting results for ■^^^ and fif 
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Figure 4: fig vs. T at = = 0.03 (left) and vs. (right) 
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We present sf^'^^ data as a function of Euclidean time in the lefthand side of Fig. Q We plot 
Bg^''^' as a function of M|- in the righthand side of Fig. Q We fit the data to the form of Eq. ( ^TTl ) 
and the preliminary results are summarized in Table ^. 

In order to perform a data analysis with higher precision, we need to know the chiral behavior 

(3/2) (3/2) 

of and By . Hence, we plan to calculate this using staggered chiral perturbation theory. 
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